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We predict a new class of excitations propagating along a Josephson vortex in 
two-dimensional Josephson junctions. These excitations are associated with the dis- 
tortion of a Josephson vortex line and have an analogy with shear waves in solid 
mechanics. Their shapes can have an arbitrary profile, which is retained when prop- 
agating. We derive a universal analytical expression for the energy of arbitrary shape 
excitations, investigate their influence on the dynamics of a vortex line, and discuss 
conditions where such excitations can be created. Finally, we show that such ex- 
citations play the role of a clock for a relativistically-moving Josephson vortex and 
suggest an experiment to measure a time dilation effect analogous to that in special 
relativity. 

PACS numbers: 05.45.Yv, 03.75.Lm, 74.50.+r 

Introduction. — The sine-Gordon (SG) model describes a variety of physical systems, 
including ferromagnets, information transport in microtubules, dislocations in crystals, non- 
linear optics, Josephson junctions and high temperature superconductors (see, e.g.,- 
The (1+1)D SG model is fully integrable and possesses well known exact solutions in the 
form of solitons and breathers^. Depending on the physical realization, they describe dif- 



ferent objects, for instance, magnetic domain walls in ferromagnets or vortices (fiuxons) in 
Josephson junctions. There exist mathematical techniques that allow constructing solutions 
of the (1+1)D SG, such as the Backlund transformation, the Hirota and inverse scattering 
methods (see, e.g., -^). Although real systems described by the SG model may contain 
additional terms associated with damping or external conditions, such terms can be taken 
into account as perturbations to the exact solutions^. 

To describe realistic systems, such as two-dimensional Josephson junctions, one should 
depart from the well studied (1+1)D model and deal with the (2+l)D SG equation, which 
is no longer fully integrable. The extra dimension opens new avenues for the existence of 
phenomena which are absent in the ID case (see, e.g.,-^ 1 ^). Here we show that there 
exists a class of exact solutions of the (2+l)D SG equation which describes the propagation 
of distortions of an arbitrary shape along a Josephson vortex line. The property of these 
excitations such as transmitting pulses of electromagnetic radiation of arbitrary shape along 
a Josephson vortex can be useful for the transmission of information in various Josephson 
devices (see, e.g., theoretical and experimental works^ 1 ^). 

Experiments^ have shown that Josephson vortices can exhibit Lorentz contraction. So 
far, another relativistic effect^, time dilation, has been considered difficult to detect in 
Josephson systems. In order to observe time dilation, an internal degree of freedom acting 
as a clock is needed. In (1+1 )D such clock is absent because solitons of the (1+1 )D SG 
equation do not possess internal oscillation modes. Recently, time dilation was proposed^ 
in junctions obeying the double sine-Gordon equation. There, an internal degree of freedom 
of a vortex, due to a special feature of the double SG model was proposed as a clock for a 
moving vortex. Here we show that there is no need to deviate from the original pure sine- 
Gordon nonlinearity because an additional degree of freedom due to an extra dimension, 
can be used to realize a clock. 

Model. — Consider a 2D Josephson junction described by the (2+l)D SG equation, 

<ftt ~ <fixx ~ Vyy + sin = (1) 

where ip is the superconducting phase difference across the Josephson junction. Here, the 
coordinates x and y are normalized by the Josephson penetration length Aj, and the time t 
is normalized by the inverse Josephson plasma frequency In the case of a Josephson 

junction, the soliton solutions of Eq. ([1]) describe Josephson vortices or fluxons^ 1 ^. 
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A more general class of solutions of the (2+l)D SG equation has the following form, 

ip(x, y,t) = 4 arctanexp [y — f(x ± t)] (2) 

which satisfies Eq. ([T]) exactly with an arbitrary real-valued twice-differentiable function 
/(£). The excitations, described by /, are associated with a local shift of a Josephson 
vortex in the transverse direction and have analogies with elastic shear waves in, e.g., solid 
mechanics and seismology^. Such distortions propagate along the stationary vortex line 
with the fixed speed \u\ = 1 in units of the Swihart velocity c = Xju p . Hereafter, we will 
refer to these excitations as shape waves. 

Because of the arbitrariness of /, Eq. (j2J) describes a variety of excitations of vari- 
ous shapes. Choosing / localized in a finite area, e.g., f = A/ cosh(x — t), Eq. (j2J) 
describes an excitation, localized along x that keeps its shape when propagating, i.e., a 
solitary wave^. For each solitary wave of this type, there exists an anti-partner with an 
/ of opposite sign in Eq.([2l). For solitary waves to be solitons, there is an additional 
important criterion: restoring their shapes after they collide. Consider a trial function 
(f(x, y,t) — 4 arctanexp [y — f(x + 1) ± f(x — t)] that, when t — > — oo, describes the propa- 
gation of two solitary shape waves toward each other (minus sign) or a solitary wave and its 
anti-partner (plus sign). One can see that (T2J) can only approximately satisfy Eq. ([1]) when 
\f'(x + t)f'(x — t)\ <C 1 for all values of x and t. This suggests that, in general, the condition 
for restoring the shapes may not be satisfied. Indeed, our numerical solutions of the (2+l)D 
SG equation show that the waves keep their shapes when propagating, but their collision is 
destructive. The collision of two waves of large amplitudes may lead to their mutual anni- 
hilation, dissipating energy away from the vortex, see Fig. [TJ However, smaller-amplitude 
waves behave similarly to solitons, dissipating very small energy and keeping their original 
shapes after collision. 

Since the equation ([1]) is Lorentz- invariant, one can obtain other solutions performing 
Lorentz transformations on Eq. ([2]). A Lorentz transformation along the x-axis leads to a 
Lorentz contraction of the excitation, but does not lead to any solution outside the class 
of Eq. (j2J), only rescaling the arbitrary function /(£). The Lorentz transformation along y 
turns out to be far more interesting, leading to a class of solutions more general than Eq. (j2J). 
After such transformation, 

<f(x, y,t) — 4 arctanexp 
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FIG. 1: (Color online) Collision and annihilation of two large-amplitude solitary shape waves 
propagating along a Josephson vortex line, Eq. ([2]) with / = 8/ cosh[(x — xo±t)/2]. The color scale 
represents the superconducting phase difference ip. The figure shows ip(x, y, t) snapshots obtained 
by numerically solving the (2+l)D SG equation on a region of size 60 x 40 Aj with Neumann 
boundary conditions: n • V</? = 0. Initially, the two solitary waves, moving at speed \u\ = 1 toward 
each other, keep their shape while propagating, but their collision at t ~ 15 leads to their mutual 
annihilation and radiation of energy in all directions. Contrary to the large-amplitude waves shown 
here, small-amplitude shape waves satisfying \f'(x + t)f'(x -t)| < 1, interact weakly during the 
collision and behave similarly to solitons. 

which describes a Josephson vortex line moving with —1 < v < 1 and a propagating shape 
excitation described by /(£). 

We will show that, apart from the Lorentz contraction^, a propagating Josephson vor- 
tex exhibits time dilation which influences the dynamics of its shape excitations, Eq. Q. 
Recently, the Lorentz time dilation of a bound half-fluxon pair has been numerically stud- 
iecU^ for a long Josephson junction with a ferromagnetic insulator. The use of the internal 
degree of freedom of a bound half-fluxon pair was proposed there to test the time dilation. 
It was found that internal oscillations of the bound pair can play the role of a clock^. A 
Lorentz reduction of their frequency was numerically observed^ when the bound half-fluxon 
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pair was moving with relativistic velocities. Here we demonstrate that, to observe the time 
dilation effect there is no need to use the double sine-Gordon equation^ but the additional 
degree of freedom of the (2+l)D SG equation can be used to realize the internal clock. 
Indeed, the shape excitation moving along a fluxon provides such a clock. In other words, 
the position of the excitation on a fluxon acts as a "minute-hand". From Eq. ([3]) the po- 
sition of a vortex changes with time as yo(t) = vt, while the position of the excitation is 
%o(t) — -f[t — v yo{t))/{l — v 2 ) 1 ! 2 = =Ft (1 — v 2 ') 1 ! 2 . Thus, the excitation moves along x with 
speed \u\ = \xq\ = (1 — v 2 ) 1 ^ 2 < 1. In other words, all dynamic processes related to a 
relativistic moving vortex are slowed down by a time dilation factor (1 — v 2 ) 1 ^ 2 , as in special 
relativity, compared to the stationary vortex. Our numerical calculations confirm that the 
collision of two solitary waves is also slowed down, compared to the collision at the rest 
frame as in Fig. [TJ 

Experimentally, shape waves can be excited when a fluxon line interacts with inhomo- 
geneities introduced on its path. When a vortex line passes through a shaped region of a 
Josephson junction, shape distortions of a vortex line are excited by inhomogeneities of the 
boundary. The shape of the excited waves is determined by the shape of the inhomogeneity 
as well as the vortex velocity. Below we present our analytical and numerical results for 
the case of a vortex line interacting with a barrier formed by a Josephson junction with 
a locally-increased width, as shown in Fig. [2^,. Another intriguing possibility to generate 
the shape excitation is to use a T junction to create a cloning barrier- on the path of a 
propagating fluxon line. Other types of inhomogeneities, besides the geometrical ones, such 
as microshorts or microresistors, realized by variation of the critical current (see, e.g.,-) can 
initiate shape excitations. One more possibility is using injectors of electric current con- 
nected locally to a Josephson junction^ where the fluxon line is propagating. This allows to 
conveniently control the shape of the excitations by varying the electric current of injectors. 

In order to detect the propagation of a shape wave we propose the following scheme which 
involves the use of a Y-shaped cloning junction^. Such junction can split a vortex into two 
parts, which afterward propagate independently from each other. For this purpose, a Y 
junction with a small span angle is preferable, in order to minimize the disturbance on the 
vortex line. When a vortex with a propagating shape wave is split by such a junction, the 
shape excitation is directed to one of the branches of the Y junction, see Fig. [2j A shape 
wave can be generated from a boundary by locally increasing the width of the Josephson 
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junction, as shown in the inset of Fig. [2^. Because of the associated energy barrier, one end 
of the fluxon line is delayed with respect to the other one, which results in a kink profile, 
as shown in Fig. [2J The generated kink starts to propagate along the vortex line away from 
the boundary, carrying information about the induced delay to the rest of the vortex. 

Consider now an arbitrary shape excitation with function / satisfying f(—W) — > — A//2, 
f(W) — > A//2 for W — * oo. The total energy of this shape excitation is the difference of the 
energy E[tp] of the solution Eq. (j3J) and the energy E[ip]f=o of the vortex without excitation: 
AE — E[ip] — E[ip]f=Q. We obtain the following analytic expression for the total energy: 



which is exact in the limit W — > oo and valid for an arbitrary excitation /(£) of height Af. As 
seen from Eq. (j4j), the energy AE can take negative values if / is a kink (Af ^ 0) propagating 
along a moving vortex line. Indeed, this is in agreement with our numerical simulations on 
Fig. [2] which indicate that the increased width AW of the Josephson junction is compensated 
by the energy gained due to the propagating kink with negative AE. Neglecting radiation 
losses, we now write the conservation of energy relation for a moving vortex line: AE + 
8 AW (1 — v 2 )^ 1 ^ 2 = 0. Using Eq. (j3J), the minimal size of the kink can be estimated: 
Af > Af m i n (AW, v) = AW /v. Note that the apparent height Ay of a kink propagating 
along a moving vortex line is Lorentz contracted: Ay = Af a/1 v 2 . To describe the kinks 
in Fig. [2] we use the approximation /(£) = |A/tanh (2v£/Af). Substituting it in Eq. (TjJ we 
obtain AE = — 8 Af v (1 — v 2 ) _1 ' 2 /3. From energy conservation, we calculate the apparent 
height of the kink to be Ay = 3 yl — v 2 AW/ v. Substituting the value AW = 2, we obtain 
the following results for the apparent heights of the kinks: Ay ~ 4.5 for v = 0.8 and Ay ~ 2.9 
for v = 0.9, which agree with the numerical results shown in Figs. [2h and [2b, respectively. 

After passing the Y junction, each half of the initial vortex line is directed to a Josephson 
junction of width W. The presence of the shape kink slows down the average speed of a 
vortex line, 



Thus, the presence of a shape wave can be directly detected by measuring voltage on the 
Josephson junction, proportional to (v). A convenient system to study shape excitations 
would be a Y junction introduced into a ring, forming a a-shaped Josephson junction which 
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FIG. 2: Detection of a vortex shape excitation and measurement of its speed in a Y junction. 
The figure shows snapshots of the superconducting phase gradient |Vy| obtained by numerically 
solving the (2+l)D SG equation in a Y junction geometry of size 40 x 80 Aj with Neumann 
boundary conditions: n • V</? = 0. At time t=0, a fluxon line moving with velocity v = 0.8 (a) or 
v = 0.9 (b) approaches a geometrical energy barrier (shown on the left inset). When passing the 
barrier, a shape excitation (zoomed on the right inset) is induced on the left end of a fluxon line 
and propagating to the right along the vortex. The propagating fluxon is split by the Y junction 
in two parts, one containing the shape excitation, and the part without it. Depending on velocity 
v of the vortex line, the shape excitation can be directed either to the right branch, if the velocity 
v is low (a), or, to the left branch, if v is higher (b). This is due to the time dilation reducting 
the speed of the shape excitation. Measuring the critical velocity v c of the vortex line when such 
a transition occurs, the speed of a shape excitation can be found as u = v c W/D. The presence of 
a shape excitation reduces the average velocity (v) of a vortex line, for instance, in (b) the vortex 
in the left branch at t = 80 is noticeably delayed from the vortex in the right branch of the Y 
junction. 
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has been recently realized experimentally with high-T c in the context of generation of Tera- 
hertz waves^i. 

The effect of time dilation is to slow down the shape excitation when the vortex line 
reaches relativistic velocities. If the velocity of the shape excitation is smaller than some 
critical value, it will be directed to the left branch of the Y junction, as on the Fig. [2b, 
while the vortex propagating in the right branch will remain undisturbed. On the other 
hand, if the shape excitation propagates faster, it will be directed to the right branch, 
Fig. [2ji. Because of the change in the average velocity of a vortex line Eq. (151) , a voltage 
step, corresponding to a shape wave changing the branch of the Y junction, will be exhibited 
on the I-V characteristics. Finding the velocity of a fluxon line v c from the position of the 
voltage step, speed of the shape excitation can be calculated as u = v c W/D, where D is 
a distance between the barrier and the Y junction, see Fig. [2j In this way, the relativistic 
factor (1 — -U 2 ) 1 / 2 for the delay of the shape kink propagation can be tested. 

In an underdamped Josephson junction, a fluxon can be accelerated to velocities close 
to the Swihart velocity, typically of the order 10 7 m/s 2 for Nb Josephson junctions. The 
velocity of a fluxon depends both on damping and an applied current^. For the damping 
parameter a = 0.01 and current 7 = 0.03 in normalized units^, it is possible to accelerate 
a vortex to the velocity v ~ 0.92, so that the relativistic time dilation factor would be 
(1 — v 2 ) 1 / 2 ~ 0.4. For instance, for typical Nb/Al— A10 x /Nb Josephson junctions with the 
plasma frequency u p /2ir ~ 50 GHz, the propagation time of the shape excitation between 
the two ends of a vortex line in a configuration as in Fig. [2] would be approximately 0.4 ns 
for a stationary vortex, while for a vortex moving with velocity v = 0.92 this value would 
be relativistically delayed by 0.6 ns. 

Conclusion. — We have predicted the existence of excitations of arbitrary shapes propa- 
gating along stationary or moving Josephson vortices in 2D Josephson junctions, and pre- 
sented the exact solutions of the (2+l)D SG equation. Remarkably, at some conditions, a 
shape wave can carry negative energy, i.e., a vortex with a shape excitation can have smaller 
energy than the same vortex without it. Based on our analysis, we suggest how to test a 
time dilation effect analogous to that in special relativity: using the shape excitation as a 
"minute-hand" measuring the proper time of the vortex's coordinate frame. 
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